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Inclusive 7r° spectra at high transverse momentum in d-Au collisions at 
RHIC 
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Preliminary results on inclusive neutral pion production in d-Au collisions at -^/snn = 
200 GeV in the pseudo-rapidity range < rj < 1 are presented. The measurement is 
performed using the STAR Barrel Electromagnetic Calorimeter (BEMC). In this paper, 
the analysis of the first BEMC hadron measurement is described and the results are 
compared with earlier RHIC findings. The 7r° invariant differential cross sections show 
good agreement with next-to-leading order (NLO) perturbative QCD calculations. 

1. Introduction 

A suppression of high hadron production relative to a simple scaling from p-p has 
been observed in central Au-Au collisions at RHIC [U It was also found that jet- 
like correlations opposite to trigger jets are suppressed and that the elliptic anisotropy in 
hadron emission persists out to very high pr [El El- In contrast, no suppression effects 
were seen in d-Au collisions, which provide an important control experiment for the effects 
of cold nuclear matter. This has led to the conclusion that the observations are due to 
the high density final state in Au-Au. The most probable explanation to date is parton 
energy loss from induced gluon radiation (jet quenching). For a recent review see Ref. 
To quantitatively understand the existing modification from cold nuclear matter, precise 
measurements of identified hadrons at high px in d-Au are required. The STAR EMC 
allows high transverse momentum measurements of ir°, 7] mesons and direct photons and 
may contribute to the identification of p mesons. In this paper we present preliminary 
results of neutral pion production in d-Au collisions. 



2. Experimental setup 

The results were obtained using two components of the STAR detector system [ E], 
namely the Time Projection Chamber (TPC) and the Barrel Electromagnetic Calorimeter 
(BEMC). The TPC is situated in a 0.5 Tesla solenoidal magnetic field and provides a pre- 
cise measurement of the charged particle trajectories. The BEMC [ 7 is a lead-scintillator 
sampling calorimeter with a of depth of 21 radiation lengths (Xq) and an inner radius of 
220 cm, divided into towers of granularity (A77, A<p) = (0.05, 0.05). Two layers of gaseous 
shower maximum detectors (SMD), located at a depth of 5 Xq, measure the EM shower 
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Figure 1. Invariant mass distribution of neutral pion candidates in minimum bias d-Au reactions 
at \J snn = 200 GeV (left). The solid line is a Gaussian plus second order polynomial fit and the 
lower histogram is obtained from event mixing. The middle (right) figure is the background subtracted 
spectrum using the polynomial fit (event mixing). The relative mass resolution is 20%. 

shape with high resolution (A77, A0) = (0.007,0.007). In this analysis a partial imple- 
mentation of the BEMC was used, consisting of 2400 towers covering < rj < 1 and 
full azimuth (the deuteron beam has positive rapidity). When complete, the BEMC will 
have a coverage of — 1 < 77 < 1. Beam test results [|H] provide the absolute energy cali- 
bration, whereas the relative calibration is obtained from the peak position of minimum 
ionizing particles (mostly charged hadrons) on a tower-by-tower basis [IHj- Moreover, an 
overall gain calibration was obtained by comparing the momentum of electrons identified 
in the TPC with the energy deposited in the BEMC. The achieved energy resolution is 
5E/E w 16%/ VE. 

In the minimum bias trigger a neutron signal was required in the Zero Degree Calorimeter 
(ZDC) in the Au beam direction resulting in an acceptance of (95 ±3)% of the d-Au 
hadronic cross section. To enhance the high range, two high tower triggers (HT1 and 
HT2) were used with an energy threshold of 2.5 GeV and 5 GeV for the highest EMC 
cluster energy, respectively. The tower occupancy is 1-5% for d-Au events, and the high 
tower trigger efficiency is nearly 100%. 

3. Data analysis 

After event quality cuts, 10M d-Au events taken in the year 2003 were used for the 
analysis presented in this paper. Neutral pions were reconstructed in the decay chan- 
nel 7T° — > 77 (branching fraction 98.8%) by calculating the invariant mass of all pairs 
of clusters in the calorimeter. Only those clusters were selected which do not have a 
TPC track pointing to them. Furthermore, a cut on the energy asymmetry was imposed 
\Ei — E2\/{E\ + E2) < 0.5. At present, the full calibration of the calorimeter and the 
finding of noisy and dead towers are in progress. To perform a neutral pion analysis at 
this stage a sub-sample of good towers was selected. The quality of the individual towers 
was checked using the ir° invariant mass distribution. The towers are assigned by the 
decay photons with the highest energy. Only those towers were used which have a well 
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Figure 2. Left: Inclusive pt distribution for neutral pions in d-Au collisions at < r\ < 1 (solid 
circles). The error bars (shaded bands) represent the statistical (systematic) uncertainties. Previous 
STAR measurements on the charged hadron cross section [E3] are shown by the open triangles. The full 
triangles show 7r° results from PHENIX [ 11 . Right: The invariant differential cross section (full circles) 
compared to NLO pQCD calculations (solid line). The factorization scale uncertainty is illustrated by 
dashed line. The lower plots show the ratios of the data to a fit of the PHENIX results (left) and to the 
theory curve (right). The dashed lines indicate the systematic errors derived from the data. The dotted 
curves (right) were obtained using different factorization scales. 



defined 7r° signal above the combinatorial background. A Gaussian fit to the tt° signal has 
to have a relative mass and width error of 30% and 50%, respectively. By this method 
approximately one third of all the towers was used for the present analysis. The tt° peak 
position of all good towers was used to perform an additional tower gain correction (7% on 
average). In Fig. ^ the resulting invariant mass distribution of cluster pairs is shown. A 
clear signal is observed with a RMS width of 28 MeV/c 2 . The combinatorial background 
from random pairs was estimated by two different methods. First, a second-order poly- 
nomial was fitted to the invariant mass distribution outside the peak region (solid curve 
in the left plot of Fig. Q). The second method describes the background using the event 
mixing technique. The event mixing distribution, which is shown in the lower histogram 
in the left plot of Fig. Q describes the background reasonably well in the range 0.8-5 
GeV/c 2 (not shown). At lower invariant mass the excess can be attributed to the tail 
of the 7T° peak and a contribution from the rj signal (m v = 547.3 MeV/c 2 ). The peak 
observed at m < 0.05 GeV/c 2 stems from cluster splittings in the EMC towers. 
The background subtracted spectra are shown in the middle and right-hand plot of Fig. ^ 
The yields per event obtained from both subtraction methods were extracted in pt bins 
(width of 0.5 GeV for minimum bias and 1 GeV for high tower triggered events) by 
integrating the background subtracted mass distribution in a range ±3cr around the ir° 
peak. The mean values were used for further analysis and the difference contributes 10- 
15% to the systematic uncertainties. The signal-to-background ratio increases from 1 at 
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about pr = 1 GeV/c to 8 at 4 GeV/c. 

Corrections for reconstruction losses (quality cuts) and detector efficiencies were calcu- 
lated with Monte-Carlo simulations using the STAR detector geometry and reconstruction 
software. A correction for the unmeasured trigger fraction, which is expected to be a few 
percent, is not applied. Losses due to cluster density effects and contributions from weak 
decays of K° mesons are not corrected yet, but are expected to be small. The high tower 
trigger spectra are normalized using pre-scale factors obtained from the ratios of the 
BEMC cluster transverse energy distributions in the overlap region. The overall system- 
atic errors related to efficiency, yield extraction, pre-scale factors, and energy calibration 
are estimated to be 30% for low and 50% for high transverse momenta. 

4. Results and discussions 

The inclusive pt distribution for neutral pions is shown in the left panel of Fig. El The 
yields up to 6 GeV/c are from minimum bias events while above 6 (9.5) GeV/c they 
are from HT1 (HT2) triggered events. For the different trigger samples the yields have 
an overlap of one point in the Pt spectrum and agree within errors. It is seen from 
Fig. |2 that 7T° mesons are presently measured up to pr = 16 GeV/c. The pr spectrum is 
compared with previous STAR measurements on charged hadron cross section [EI] and 
with PHENIX results on neutral pion production [ ■ There is a reasonable agreement 
within 10-20%. 

The right panel of Fig. |21 shows the invariant differential cross sections obtained from 
the product of the yields and the hadronic cross section in d-Au collisions [ (TUJ . The 
normalization uncertainty is 10%. The results are compared to next-to-leading order 
(NLO) pQCD predictions [ H2] which are calculated using the CTEQ6M set of nucleon 
parton distribution functions [ [QJj and the nuclear parton densities in the gold nucleus 
from [E] • I n this calculation the factorization scale was identified with pt and is varied by 
a factor two to estimate the scale uncertainties (see Fig. EJ). The fragmentation functions 
are taken from The Cronin effect is not included in the calculations. The data, 

within errors, are consistent with the calculation up to p? = 15 GeV/c. 
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